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INTRODUCTION 


Traditionnally high purity silicon for solar and electronic applications 
has been produced in Siemens decomposers by hydrogen reduction of chlorosilanes 
[1], The yield in the Siemens process is poor. The new cheap routes to silane 
offer an attractive alternative for silicon production via silane pyrolysis. 

The thermodynamic yield of this process is 100%. However, higher feed concen- 
trations of silane (a few percent) in Siemens decomposers lead to formation of 
fines via homogeneous nucleation. In order to realize the potential of silane 
pyrolysis the formation of fines must be prevented or at least kept below a 
certain acceptable minimum level. Fluidized beds offer an attractive possibil- 
ity, and extensive experimental studies have been conducted [2,3]. Modeling 
efforts have lagged far behind the experimental activities. 

The objective of this paper is to develop a mathematical model for fluid- 
ized bed pyrolysis of silane that relates production rate and product silicon 
properties (such as size, size distribution, presence and absence of fines) with 
fluidized-bed size and operating parameters (such as wall temperature, feed 
concentration, gas flow rate, seed size, etc.) and with bed grid design. Upon 
model verification it is desired to expand the model to account for product 
morphology, i.e. porous versus nonporous particles, etc. While fluidized-bed 
models for catalytic processes are abundant (e.g. see reviews by Grace or Yates 
[4,5,6]) and models for gas-solid reactions with changing solids have been well 
established by Kunii and Levenspiel [7,8], a comprehensive model needed for 
simultaneous chemical vapor deposition (CVD) and nucleation reactions has not 
been reported. 


MODEL DEVELOPMENT 

A suitable mathematical model for silane pyrolysis in fluidized-bed reac- 
tors should consider various reaction pathways, the problem of "smoke" (fines) 
formation, the suppression of "smoke" formation and its capture by large parti- 
cles. We approach these problems on two levels. First, we attempt to identify 
a plausible description of the key chemical and physical rate determining steps 
in reaction pathways from silane to silicon. Then we address the question of 
flow and gas-solid contacting in fluidized-beds, and develop a model for the 
reactor. We treat here only batch growth of solids because data are readily 
available only for this model of operation. 

Reaction Pathways from Silane to Solid Silicon 

The mechanism of silane pyrolysis is not completely understood [1]. Active 
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work is in progress in addressing this question [9*10] . We do not know for 
certain which intermediate species is the first to nucleate. Based on the 
current understanding of the system we arrive at the picture presented in 
Figure 1 for the various pathways from silane to silicon. We assign the rate 
of the rate limiting step to each pathway. In doing so we have assumed that 
silane can decompose by two independent pathways. One is the homogeneous de- 
composition (pathway 3) into a gaseous precursor that can nucleate a new solid 
phase of silicon. We use the form reported by Hogness et. al. [11] with the 
variations suggested by either Purnell and Walsh [12] or O'Neal and Ring [13] 
to describe it. The other route for silane decomposition is the heterogeneous 
chemical vapor decompostion (CVD) of silane on the existing silicon seed parti- 
cles (pathway 1) or on the formed nuclei (pathway 2) . For this rate we use the 
first-order form reported by Iya et. al. [14]. It can readily be shown that in 
fluidized beds, due to rapid particle-gas transfer, CVD would control the 
growth rate of large particles rather than mass transfer. Since the molecular 
bombardment rate of small particles (fines) is larger than the CVD rate at the 
temperature of interest (500 to 750°C) then CVD also controls the silane loss 
to fines (pathway 2). 

The intermediate which leads to the nucleation of the new phase is assumed 
to be the silicon vapor. The concentration of this intermediate formed by 
homogeneous pathway 3 is always very small. By pathway 4 nucleation of criti- 
cal size nuclei, r*, occurs whenever supersaturation S = 1 is exceeded. The 
concentration of silicon vapor can be suppressed by diffusion and condensation 
on large particles (pathway 6) and by molecular bombardment of fines (pathway 
5) . We assume here that nucleation occurs by the homogeneous nucleation theory 
[15]. The molecular bombardment rate of small particles (pathway 5) is calcu- 
lated by the classical expression of kinetic theory [16] while the diffusion 
rate to large particles (pathway 6) is readily obtained from film theory of 
mass transfer [16]. This concludes the description of those pathways in Figure 
1 that include gaseous and solid' species . The rate forms used are indicated on 
the figure. 


Population Balance for Fines 


In order to describe the size of fines and the ability of large particles 
to scavenge them, we need a population balance for fines. We define n v (v)dv to 
be the number of fines having volumes in the range v to v + dv per unit volume 
of the gas. The population balance, based on a well mixed gas volume, yields 
the following equation: 


V^ n v (v) + h [R G (V) n v (v)] = \ ; v* 3 ( v ~v,v) n v (v-v) n v (v)dv 

- / B(v,v) n (v) n (v)dv - a(v,v T ) n (v) + r u . T N. 6 (v-v*) (1) 

V V L V riJN A 


The LHS terms of Eq. (1) are the elutriation rate and growth rate into size 
range by CVD of silane and molecular bombardment of Si vapor. The RHS terms are 
growth rate into size range by coagulation of fines, rate of loss in size range 
by coagulation, rate of loss in size range by scavenging by large particles 
(seeds) and the nucleation rate of the critical size. . The pseudo -state approxi- 
mation is used here since we are interested in observing the system at the time 
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Coagulation and coalescence of fines— eq. (4) 

Scavenging by seed particles on fines - eq. (10) 

Attrition of large particles (not Included) 


Figure 1. Schematic of Various Pathways for Conversion of 
Silane to Silicon. 
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scale for growth of large particles which is very long compared to other time scales 

It should be kept in mind that we are not interested in a precise descrip- 
tion of the size distribution of fines. We mainly need to establish if their 
formation can be suppressed and if not, what their mean size might be. We, 
therefore, solve the population balance approximately by using the method of 
moments. We define: 

M» = /°°. n (v) dv = (total number of particles per unit (2) 

0 v v volume of gas) 

00 

M = / ^ vn^(v)dv = (total volume of particles per unit (3) 

V volume of gas) 


The number average mean volume of fines is then obtained by the ratio of the 
first and zeroth moment v = Mp/Mg- In order to reduce eq. (1) to a set of equa- 
tions for the moments we need to make a set of assumptions. These are: 

i) The coagulation coefficient is assumed to be a constant which only depends 
on the average size of fines. This approach is often used in the continuum 
regime where the problem is thus reduced to a Smoluchowski type equation [16,17]. 
We use the same approach for the free molecular regime and take the following 
value of the coagulation coefficient Bo 

continuum region (4a) 

free molecular regime (4b) 


8kT 
"0 l3y 


5/ 2 ,3, 1/6 ( 6kTv -1/6 
^4tt 'p ' v 


ii) The scavenging coefficient is assumed to be dependent only on the mean 
size (diameter) of the large seed, particles and of the fines. Following the 
work of Doganoglu et. al. [18] and Peters et. al. [19] we represent the scaveng- 
ing coefficient of fines by large particles in a fluidized bed by: 


a(v L'V " E “mf 


3 

23 PL E mf 


(5) 


where £ m f, u m f are bed voidage and superficial gas velocity at minimum fluidiza- 
tion, a is the scavenging coefficient, dpp is the mean diameter of large parti- 
cles (seeds) and E is the single large particle collection efficiency. It is 
well established [16,17] that single body collection efficiency can be approxi- 
mately represented as a sum of efficiency for impaction, interception, diffusion, 
and diffusion with interception. The diffusional mechanism dominates under the 
conditions prevailing in a fluidized bed reactor and 

E = 2 Pe“ 2/3 (6) 

where Pe = dpp u m f/D and D = kT/3iry dpp. The mean particle diameter for small 
and large particles is estimated. 

iii) The growth rate of fines, which is due to the combined effect of CVD 
growth from silane and molecular bombardment of gaseous silicon, is assumed to 
be given as a function of the average size of fines, i.e . 
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where o = (M sl /<f> g P s± ) (4 ) 1 / 3 3 2 / 3 (r HT + r^) . 


With the above three assumptions the population balance, eq. (1), is 
reduced to the following two equations for the moments: 


2 e 0 M 0 (0t + Vj* M 0 " r HN N A ” 0 


( 8 ) 


a m, 


1/3 M 2/3 

M, 


0 


- (a + -) M 1 + r HN N A v^ = 0 


(9) 


The scavenging rate of fines by large particles is now given by: 
Si 
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with the scavenging coefficient a given by eq. (5), 


This completes the description of various pathways of Figure 1, namely 
pathway 7, by which a population of fines with a certain average size is estab- 
lished, and of pathway 8 for scavenging (filtering) of fines by large particles. 

Backmixed Reactor (CSTR) Model 

The hydrodynamics of a fluidized bed and gas-particle contacting are com- 
plex and not entirely understood. It is useful, however, to develop models that 
characterize the limiting behavior of the system which can only be approached in 
practice. Here we deal with the competitive homogeneous and heterogeneous reac- 
tion. Intimate gas-solid contacting, high solids to gas ratio and no gas by- 
passing will favor the heterogeneous route. Homogeneous nucleation is to be 
prevented compared to CVD growth and diffusion of Si vapor. The former can be 
regarded as a reaction of high order and, hence, will be suppressed the most by 
complete micromixing such as found in an ideal CSTR. Therefore, both suppres- 
sion of homogeneous decomposition and of homogeneous nucleation will be favored 
in a CSTR, i.e . in an ideally backmixed reactor. This situation can be ap- 
proached in fluidized beds when bubble formation is suppressed while good solids 
mixing is maintained. 

We formulated the CSTR model based on the following assumptions: i) no 

wall deposition, i.e . negligible wall to particles area, ii) no temperature 
gradients between the gas and particles, iii) uniform composition and tempera- 
ture in the reactor. In addition, for the results presented in this paper, we 
assumed that all seed particles have the same initial size and grow at the same 
rate . 


The equations to be solved are: 
i) The mass balance on silane: 
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ii) The mass balance on silicon vapor: 
0 


q C + Vr 

e Si,e HD 


Vr HN + A TL r DL + A TF r DF 


( 12 ) 


iii) The balance on fines given by eqs. (8-9): 
iv) The energy balance which can be shown to readily simplify to the follow- 
ing equation for temperature: 


T = 
e 


h A/h A 
DP w w 

1 + h D A D /h A ~D 


T„ + 


1 + h^A^/h A 
D D w w 


w 


(13) 


w w D 

v) The growth rate of large seed particles 

dR L M Si ^ r HT + r DL^ 


m 


dt 


Si 


<J> 


+ 


Sea 


(14) 


4 1,R 1 *L P S1 


The total surface area of large and small particles ApL, App, respectively, are 
estimated from their average size. 

A computer program was developed for solution of eqs. (11-14) and eqs. (8- 
9). In addition to the quantities already discussed many parameters of the 
fluidized bed such as bed height at minimum fluidization, the elutriation con- 
stant, heat transfer coefficients, etc. are estimated from the literature [7,20, 
21]. It should be noted that the distributor must be cooled to keep its temp- 
erature below 350° to prevent silane CVD in the nozzles. When the distributor 
and wall temperatures are known bed temperature is given by eq. (13). When bed 
temperature is known eq. (13) is bypassed in the program. 

The required input variables for the CSTR model are bed diameter, initial 
weight of solids, initial size of solids, gas flow rate, inlet gas temperature, 
inlet gas composition and pressure. Equations (8-9), (11) and (12) are solved 
for the initial seed size and the RHS of eq. (14) is then evaluated. The pro- 
cess is repeated for various selected values of seed size. This generates the 
set of pairs of values of Rp vs dRp/dt which is numerically integrated to estab- 
lish the relationship between time and particle (seed) radius. It should be 
noted that the CSTR as presented here has no adjustable parameters. Given the 
input quantities, the rate forms and transport properties are calculated by the 
appropriate subroutines and reactor performance is predicted, i.e . the growth 
rate of seeds, amounts of fines elutriated and silane conversion are calculated. 


The fraction of silane that ended in the form of free fines, Ff, is calcu- 
lated at each set of conditions from 


F 


f 


1 


growth rate of large particles 
total rate of silane deposition 


1 - 


A TL( r HT +r DL )M Si +m Sca 

( ' A TL r HT +A TF r HT +Vr HD' )M Si 


(15) 
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Bubbling Fluidized Bed (FBBR) Model 


A large diameter fluidized bed of silicon particles operated at three to 
ten times the minimum fluidization velocity will behave like a bubbling bed 
provided bed diameter is large enough. However, since the information regard- 
ing the kinetics of the system under study is not very refined, a sophisticated 
bubbling bed model is not necessary. Only the key features of the bubbling bed 
need to be incorporated in the model. These are the existence of three dis- 
tinct regions: the grid region, emulsion region (including the cloud) and the 
bubble region. We assume that in beds of large particles jets form in the grid 
region. We assume gas to be in plug flow in the jets, and jet penetration 
height is calculated from the Yang and Realms' [22] formulas. When jets break 
up, gas is assumed to pass through the emulsion phase at minimum fluidization 
velocity. There are no gradients between gas and particles in the emulsion 
phase which is assumed well mixed. The excess gas, u-u m f, forms bubbles which 
rise in plug flow. Average bubble diameter, average bubble volume fraction, 
bubble rise velocity, etc. are calcualted from the available correlations. The 
overall mass and heat transfer coefficients between bubbles and emulsion are 
evaluated based on the approach suggested by Kunii and Levenspiel [7]. The 
mass and heat transfer exchange coefficients between jets and emulsion are cal- 
culated as a multiple of bubble-emulsion exchange coefficients as suggested by 
Weimer and Clough [23]. This multiplicative factor fj^ defined by 


K. 

je 



^be'b 


(16) 


is the only adjustable parameter for the bubbling bed model. It is also 
assumed that no large particles are present in the jets and bubbles. 


The schematic of the bubbling bed model is shown in Figure 2. The various 
pathways described in Figure 1 can take place now in each of the three regions: 
jets, bubbles and emulsion. In addition, one must account for the interchange 
between the various regions. The equations for the mass balance on silane, 
silicon vapor and fines and for the energy balance in each of the three regions 
are lengthy and will be omitted here but are given elsewhere [24]. 


The input parameters besides those required by the CSTR model include the 
specification of the distributor, i.e . number or orifice holes, hole diameter, 
etc. The details of the computer algorithm are reported elsewhere [24], Given 
all input parameters the program calculates silane conversion, fraction of 
silane in form of fines and growth rate of large particles. Repeated calcula- 
tions at various particle sizes lead to desired time on stream-particle size 
relationship . 


RESULTS AND DISCUSSION 


Model predicted results are compared with two JPL experimental runs in 
Table 1. Both models predict reasonably well the production rate and final 
particle size. However, model discrimination is impossible because the experi- 
mental runs were of insufficient duration so that the actual final particle 
size could not be determined accurately. This also leads to problems in the 
mass balance for the experimental runs and slight discrepancies in reported 
production rates and final particle size. 
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(5) Brownian diffusion of Si Vapor 


into seed particles 


Figure 2. Schematic of the Fluidized Bubbling Bed Reactor (FBBR) 
Model for Silane Pyrolysis. 
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Table 1. Comparison of the CSTR and FBBR Model Predictions and Experimental 
Results for Two JPL Runs 


Silicon Seed Experimental Conditions 


Weight d 

Silane Feed Cone. Bed temp. Total gas 

Duration 

(Ks> (u»> 

(%) 

(°c) 

flow rate 
(moles /min) 

(min) 

10.50 227 

20 

650 


3.0 

90 

11.34 212 

80 

650 


2.5 

173 



PRODUCT COMPARISON 



Experimental 

Data 

Model Predicted 

(CSTR) 

Model Predicted (FBBR) 

Production rate 

d 

Production rate 

d 

Production 

rate d 


P 

(ym) 

(Kg/hr) 

P 

(ym) 


(ym) 

0.87 

235.5 

1.00 

237.4 

0.93 

236.6 

3.50 

241.5 

3.30 

260.3 

3.16 

257.7 


Reactor Specifications: Bed diameter 15.4 cm (6.065" I.D.); Number of orifice 

holes in distributor: 4,500. 

2 

Orifice area: 0.02 cm ; Distributor temperature: 200°C; Entering 

gas temperature: 200°C. 


In order to illustrate the dominant pathways in silane conversion to sili- 
con the modeling results for the ( two runs of Table 1 are presented schematically 
in Figure 3 for the CSTR model. It is clear that in an ideal CSTR gas-solid 
contacting is very efficient, micromixing is excellent and homogeneous nuclea- 
tion and fines formation can be effectively suppressed. For example even at 
80% SiH 4 in the feed (Figure 3b) 84% of decomposed silane (2.82 kg/h) reacts 
by CVD on the growing seed particles, 15% (0.49 kg/h) decomposes homogeneously 
and 1% (0.035 kg/h) reacts by CVD on the fines. Silane conversion is over 99%. 
The Si vapor is formed at a rate of 0.49 kg/h, but 75% of it (0.36 kg/h) is 
effectively scavenged by diffusion to and condensation on seed particles, while 
only 25% (0.12 kg/h) contributes to the mass generation of fines. The nuclea- 
tion rate is kept at a very low level of 4 x 10~5 kg/h. The nucleated fines 
gain 78% of their mass (0.12 kg/h) by molecular bombardment of silicon vapor 
and 22% (0.035 kg/h) by CVD of silane on fines. Most importantly due to 
excellent contacting 70% (0.11 kg/h) of the fines formed are scavenged by large 
particles and only 30% (0.04 kg/h) are elutriated. This means that in an ideal 
CSTR at 650°C bed temperature and at a high production rate of 3.3 kg/h only 
1.5% of silane would end in the undesirable form of fines at the reactor exit 
even at 80% SiH^ in the feed. 

Figure 4 illustrates the FBBR model predictions for the JPL run listed as 
example 1 in Table 1 at 20% SiH^ in the feed. The values computed in Figure 4 
are based on the assumption that the exchange between jets and emulsion is fifty 
times faster than between bubbles and emulsion, i.e . fj^ = 50. Homogeneous 

decomposition is favored over CVD on seeds in the jet region due to poor gas- 
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Figure 3a. Net growth rate: 1.00 (kg/h) 

Backmixed (CSTR) Reactor: 20% SiH^ in the Feed 


A TL r HT 


2.824 (kg/h) 



V I r HD 


0.488 (kg/h) 


^l/DL 


0.364 (kg/h) 


A TF r HT 


0,035 (kg/h) 


0.111 (kg/h) 



Elutriation (negligible) 

-> 


j4>(kg/h) | 3.97x10~ 5 • (kg/h) 


0.0477 (kg/h) 
F f = 1.4% 


Figure 3b. Net growth rate: 3.30 (kg/h) 

Backmixed (CSTR) Reactor: 80% SiH^ in the Feed 


Figure 3. Rates of Various Pathways for Silane Decomposition 
as Calculated by the CSTR Model for Bed Temperature 
of 650°C 
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Figure 4. Rates of Various Pathways for Silane Decomposition as 
Calculated by the FBBR Model (Run 1, T g “ 650°C, 

fjb = 50 >- 
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solid contacting in spite of the high exchange coefficients. Most of the 
homogeneous nucleation takes place in the jets (grid region) (2.6 x 10~5 kg/h) 
followed up by nucleation in the bubbles (2.5 x 10“ 5 kg/h) while little 
occurs in the emulsion phase (3.2 x 10“^ kg/h). The fines gain the most mass 
in the bubbles at a rate of 0.11 kg/h, followed by the grid region with the 
rate of 0.010 kg/h and the emulsion region with the rate of 0.0059 kg/h. 

Bubbles are the main culprit in formation of fines since 88% of the mass of 
fines is generated there. The predicted fraction of silane that is 
elutriated as silicon fines is 8.1% which exceeds the experimental value of 
3%. However, in the run of such short duration the mass balance is difficult 
to close and it is questionable whether all the fines have been detected in 
the experimental run. There are reasons to believe that actual elutriation 
rate of fines is higher than reported values. 

The comparison of CSTR and FBBR (fjb = 50) predictions for deposition 
rate with 20% SiH^ in the feed and as a function of bed temperature and gas 
flow rate is illustrated in Figure 5. At low bed temperatures and at low 
gas flow rate (i.e. close to minimum fluidization conditions) CSTR and FBBR 
predictions are close. While mathematically this is to be expected, it is 
doutbful that CSTR behavior can be achieved in practice since solids circu- 
lation would not be vigorous enough at close to minimum fluidization con- 
ditions. It is clear that the CSTR model gives an upper bound on deposition 
rate as argued earlier. At higher temperatures FBBR predictions can deviate 
significantly from CSTR behavior due to increased production of fines. This 
is to be expected for two reasons. The homogeneous decomposition has a 
higher activation energy than the CVD reaction and is favored at higher 
temperatures. At the same time bubble expansion is more drastic at higher 
temperatures and the gas bypassing problem is aggravated. 

A limited parametric sensitivity study of the model was performed. The 
changes in deposition rate, silahe conversion and formation of fines were 
determined as a function of the following quantities: (i) the kinetic form 

for homogeneous decomposition [11,12,13], (ii) the jet-emulsion exchange co- 
efficient, i.e . variation in fjb values, (iii) grid design, (iv) method of 
treatment of the population balance. It was shown [24] that the variations in 
the kinetic forms and method of treatment of the population balance for fines 
have a limited effect on production rate and fines elutriation. The model is, 
however, most sensitive to the jet-emulsion exchange coefficient. Table 2 
illustrates that an increase in the jet-emulsion exchange coefficient in the 
limit leads to CSTR behavior. Unfortunately, there is no reliable model for the 
grid region based on which the exchange coefficients could be tied initimately to 
grid design. Therefore, two different assumptions are made in order to 
estimate the effects of grid design on reactor performance. In the first 
case, it is assumed that jet emulsion exchange is governed mainly by bed 
dynamics which is primarily affected by bed diameter, height and total gas 
flow. the ratio of jet-emulsion and bubble-emulsion exchange was set a 
fjb = 50 for all grids. The effect of grid design is illustrated in Table 3. 
Clearly, silane overall conversion is hardly affected by grid design. Higher 
deposition rates (and lower formation and elutriation of fines) are obtained 
at lower jet velocities. Higher jet penetration at the same gas velocity 
also favors improved deposition rate. Yang and Keairns [22] indicate that 
the jet penetration length is proportional to d o 0-8 Uo 0.36 where d D is orifice 
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DEP. RATE (KG /HR) 



Figure 5. Comparison of CSTR and FBBR Model Predictions for 
Deposition Rate (Run 1) . 
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Table 2. Comparison of 

the FBBR Model 

at Different 

Levels of Jet- 

Emulsion 

Exchange and 

the CSTR Model 

(Conditions of 

Run 1, Table 1 

.) 



FBBR 


CSTR 


fj b = io 

fjb = 50 

fjb = ioo 

- 

Silane Conversion (%) 

99.50 

99.21 

99.15 

99.04 

Fines Elutriation (%) 

30.9 

9.1 

1.6 

0.0 

Deposition Rate (kg/h) 

0.696 

0.922 

0.986 

1.002 


Table 3. Effect of Grid Design (Constant fjb = 50) 



Grid 1 

Grid 2 

Grid 3 

Grid 4 

Grid 5 


Aj- o =0 . 002cm^ 

A ro =0. 02cm2 

A ro =0 . 02cm^ 

A ro =0. 2cm^ 

A ro =2 . Ocm^ 


N t =4, 500 

N t =450 

ft 

II 

•£> 

Ui 

N t =45 

N t = 45 

Jetting height (cm) 

0.35 

0.83 

1.97 

2.12 

2.29 

Jetting velocity 
(cm/ s) 

10.1-29.5 

101-256 

1014-1680 

101-284 

10.1-34.3 

Silane Conversion 
(%) 

99.2 

99.5 

99.7 

99.3 

99.3 

Fines Elutriation 
(%) 

8.1 

30.7 

47.8 

21.0 

0.9 

Deposition Rate 
(kg/h) 

0.92 

0.70 

0.53 

0.79 

1.00 

Nj-0 • 64A ro (cm2) 

4.4 

1.0 

0.23 

2.3 

23 


diameter and u 0 is the velocity at the orifice. At fixed total gas flow rate 
it can then be shown that the gas residence time in the jetting region is 
proportional to Nj. 0.64^2.08 = N t 0.64A. ro w here N t is the total number of 
orifice holes and A ro is the surface area of an individual orifice. Last 
row of Table 3 gives the value of this group for each grid. It is apparent that 
deposition rate correlates at least qualitatively with the residence time in 
the jetting region. 

On the other hand one can assume that the jet-emulsion exchange is 
dominated by the grid hydrodynamics and that fjfo increases proportionately 
with jet velocity. The deposition rate for grids 2, 3 and 4 which now have 
f-jb = 150, 450, 150, respectively, is increased compared to values in Table 
3, but not sufficiently to approach the superior performance of grid 5. 

However, grid 4 would yield now a deposition rate of 0.97 (kg/h) and fines 
elutriation of 3.1% which is better than the grid 1 performance. This in- 
dicates, as expected, that N t ®‘®^^A ro is not the only measure of grid 
performance. The exchange coefficient is another key variable. Its a priori 
predictions at present are not possible. Grid 5 design is also unrealistic. 

The above results indicate the importance of grid design and gas-solids 
contacting in the grid region on reactor performance. Good jet penetration, 
high gas residence time in the jetting region and excellent jet-emulsion 
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exchange are necessary for suppression of formation of fines. The current 
models for the grid region are insufficient to fully quantify reactor per- 
formance and can benefit from further improvements. It should be noted that 
the current FBBR model is capable of predicting temperature gradients in the 
grid region and that such gradients have been observed experimentally. 


CONCLUSIONS 


An ideal backmixed reactor model (CSTR) and a f luidized-bed bubbling 
reactor model (FBBR) have been developed for silane pyrolysis. Silane decom- 
position is assumed to occur via two pathways: homogeneous decomposition and 

heterogeneous CVD. Both models account for homogeneous and heterogeneous 
silane decomposition, homogeneous nucleation, coagulation and growth by 
diffusion of fines, scavenging of fines by large particles, elutriation of 
fines and CVD growth of large seed particles. At present the models do not 
account for attrition. 

The preliminary comparison of model predictions with JPL experimental 
results shows reasonable agreement. The CSTR model with no adjustable para- 
meter yields a lower bound on fines formed and upper estimate on production 
rate. The FBBR model overpredicts the formation of fines but could be 
matched to experimental data by adjusting the unknown jet-emulsion exchange 
coefficients. In the limit of low gas flow-rate and large exchange co- 
efficients the FBBR model is reduced to the CSTR model provided the jet region 
becomes negligibly small. 

The models indicate clearly that in order to suppress formation of fines 
(smoke) one must achieve good gas-solid contacting in the grid region and 
eliminate or suppress the formation of bubbles. 
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NOMENCLATURE 


a D 

A^p 

A TL 

A w 

c SiH 4 
c S,f > c S,e 


^Si,e> 

c gio 

CSi 


c Si 


d PL 


- area of the distributor excluding orifice holes, cm^ 

- total surface area of fines in the reactor, cm^ 

- total surface area of seed particles in the reactor, cm^ 

- area of vessel wall, cm^ 

- concentration of SiH 4 , mol/cm d 

- concentration of SiH 4 in the feed fas and in the backmixed 
reactor, respectively, mol/cm d 

- concentration of Si vapor, mol/cm^ 

- equilibrium Si vapor concentration, mol/cm d 

- the supersaturated concentration of silicon in the gas phase 
from homogeneous decomposition, mol/cm d 

- diameter of seed particle, cm 
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d PL> d PF 
D 

E 

f jb 
Ff 

h D 

h w 

k 

K j e 
(K b e) b 


m Sca 

M 0 

M]_ 

M Si 

n v (v,t) 


% 

Nl 

Pe 



q e 

qf 

r* 


r DF 

r DL 

r DH 

r HN 

r HT 


Rg( v ) 

S 

t 

T>T e 

t D 

u mf 

v,v F 

V L 

v* 


Vi 


- average diameter of seed particles and fines, respectively, cm 

- particle diffusion coefficient, cm 2 /s 

- single seed particle collection efficiency, Eq. (13) 

- ratio relating j etting-emulsion intercharge to bubble-emulsion 
interchange 

- fraction of silane decomposed into fines 

- heat transfer coefficient between bed and distributor, 

J/cm 2 sK 

- heat transfer coefficient between bed and surface, J/cm 2 sK 

- Boltzman constant, 1.38066 x lO -2 ^ J/molecule K 

- interchange coefficient between jets and emulsion per unit 
volume of jets, s" d 

- interchange coefficient between bubbles and emulsion phase 
based on volume of bubbles, s _d 

- mass of molecule, g/molecule 

- total amount of fines captured by seed particles, g/s 

- total number of particles per unit volume of gas, cm - 3 

- total volume of particles per unit volume of gas, cm-^/cm^ 

- molecular weight of silicon, g/mol 

- size distribution density function of fines per unit volume 

of fluid, cm +2 cm + 3 

- Avogadro's number, 6.02 x lO 2 ^ molecules/mol 

- total number of seed particles in the reactor 

- Peclet number 

- partial pressure of Si vapor, atm 

- equilibrium Si vapor pressure, atm 

- volumetric flow rate of gas in the backmixed reactor, cm 2 /s 

- inlet volumetric flow rate of gas, cm^/s 

- radius of the critical nucleus, cm 

- rate of molecular bombardment of Si vapor on fines, 
mol/cm 2 s 

- rate of molecular diffusion of Si vapor to seed particles, 
mol/cm 2 s 

- rate of homogeneous decomposition of SiH 4 > mol/cm 2 s 

- rate of homogeneous nucleation, mol/cm^s 

- rate of heterogeneous CVD of SiH^, mol/cm 2 s 

- gas constant, 8.31441 J/mol K 

- radius of seed particles, cm 

- rate of fine particle growth from heterogeneous CVD growth 
and Si molecular bombardment, cm-Vs 

- supersaturation ratio, Psi/Psi 

- time scale, s 

- temperature of the interstitial gas, K 

- temperature of distributor, K 

- temperature of vessel wall, K 

- superficial fluid velocity at minimum fluidizing conditions, 
cm/ s 

- volume of fine particles, cm 2 

- volume of seed particles, cm^ 

- critical volume of fine generated from homogeneous nucleation, 
cm 2 

- total volume of interstitial gas, cm 2 
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-e- qiq 


Greek Symbols 


a c 

a (v,v L ) 

8 (v,v) , 3o 
6 (v-v*) 

e mf 

y 

PSi 


s 


condensation coefficient 
scavenging coefficient, s - -*- 
coagulation coefficient, cm^/s 
Delta function, cm - ^ 

void fraction in a bed at minimum fluidizing conditions 
viscosity of gas, g/cm s 
density of silicon, g/cm^ 

rate constant of fine particle growth, cm/s 
specific surface free energy, J/cm 2 
sphericity of a particle 
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DISCUSSION 




LORD: Have you looked at the characteristics of an ideal grid design for 

suppressing fines? 

DUDUKOVIC: No. 

LORD: Can you give any indications of the direction which you should go, 

based on your model? 

DUDUKOVIC: Well, I am hesitant to speculate on that, because grid design was 

not part of my expertise when we started on this project. Now, after we 
have finished the study, we are convinced that it has a determining 
effect. First, we would like to crystallize our ideas as to how we should 
do an experimental study. Since it is very difficult for us to run a 
silicon system, we would use a mock system to determine the exchange 
coefficiencies and then, when we are sure that we are on the right path, 
we could go further. 

FLAGAN : Is it true that you neglected the diffusional resistance between the 

gas and the large particles that were growing? 

DUDUKOVIC: That is correct. The diffusional resistance is small compared to 

the chemical vapor deposition growth rate, and so the major resistance 
(i.e., the rate limiting step) is growth by chemical vapor deposition. 

FLAGAN: Is that true throughout the entire range of conditions that you’ve 

looked at? 

DUDUKOVIC: It seems that way. 
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